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Abstract 


In this paper we discuss a theory of normal modes of oscillation of compound liquid drops 
and the experiments performed to determine its validity. The modes are characterized by 
their frequency, the attendant displacement of fluid boundaries, and the flow pressure 
fields within the fluids. The drops consist of three fluids; a core fluid, a fluid shell 
surrounding the core, and a host fluid surrounding the shell. These fluids are assumed to 
be inviscid and incompressible, and the core and the shell to be concentric. The theory is 
obtained by linearization of the equations of fluid motion to the lowest order of nonlinear- 
ity that yields the normal modes. Numerical values of mode frequencies and the associated 
relative displacements of the fluid boundaries are presented for several specific systems, 
and the results compared with our observations. The core-centering phenomenon whereby the 
oscillations of the system tend to drive the shell and the core to be concentric was 
observed in the experiments and will be fully analyzed in a sequel. 

Introduction 


This is the first in a series of reports on the study of compound liquid drop systems. 

The systems consist of three fluids: a host fluid infinite in extent, that surrounds a 
second ^luid in the form of a shell, which in turn surrounds a third fluid that forms the 
core. In particular we will investigate the several normal modes of compound drops on the 
assumptions: (i) the density of the fluids and the interfacial tensions are arbitrary, 

(ii) the fluids are incompressible, (iii) they are inviscid, (iv) the equations of fluid 
motion are linearized to the lowest order of nonlinearity that yields the class of normal 
modes being studied, and (v) the two fluid interfaces are nearly spherical and concentric. 

In ensuing reports we relax conditions (iii) through (v) . 

The primary aim of the entire study is to gain sufficient understanding of the behavior 
of compound drops to plan and interpret experiments in the laboratories, in the weightless 
environment provided by flight on the KC-135 aircraft and the Space Shuttle. 

Aside from its interest as a fundamental study of compound drops, the work can be applied 
to the fabrication of fusion target pellets, development of container less materials pro- 
cessing techniques both terrestrial and extra-terrestrial, and development of techniques for 
liquid drop control that can be used for fundamental studies in other scientific disciplines 
such as superfluid drop dynamics. 


Equations of motion 


We consider a system of several inviscid, incompressible fluids that are in contact with 
one another at the fluid boundaries each of which is characterized by an interfacial ten- 
sion. The equations of motion fo. the system are well known and in spherical coordinates 
are given as: 



( 1 ) 

( 2 ) 

(3) 


The first equation for the velocity potential in each fluid follows from the assumption 
of incomDressibility . The second equation is also kinematic and states that the fluid 
boundary moves with the fluids; the equation of a boundary being given as r = R(0 , , t) . In 
eq. (2), f denotes the unit vector along r, the vector from the origin to R(e, 4 i), and Vg 
denotes the surface gradient operator given by 
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Equation (3) is a dynamical equation indicating that the stress across an interface 
depends solely on its interfacial tension, o. The left side of eq. (3) is the stress at a 
point on the boundary, n is the surface normal at the point, and Vg denotes the surface 
divergence operator. The right side of this equation is simply the pressure difference sup- 
ported by the boundary stress, and P denotes the pressure in the resting fluid. 

These nonlinear equations are rendered linear by assuming that the motion of the system 
produces small deviations of a boundary from the shape it had when the system was at rest. 
For a compound drop the boundaries are spherical, and the pressure difference across a 
boundary of radius f[ and interfacial tension, o, is 
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For this system, retaining only first order terms, eq. (2) can be written as 


3t ■ 3r 


0 , 


(2A) 


and eq. (3) after substituting eq. (5), becomes 

A (p II) = ^ (2 - L^) AR , (3A) 

_ 2 

where aR = R - R, and L is the surface Laplacian which results from the linearization of 
the curvature. 


Normal modes of a concentric three fluid systems 

We consider a concentric three fluid system as shown in Figure 1 , and calculate normal 
modes using eqs. (1), (2A) , and (3A) . The solution of the equations involves finding flow 
potential, , in the three fluids. Since is nonsecular solution of Laplace equation, it 
can be expressed in the following forms in each region of the system: 


:*((£, t) = ^ [A(i,mst)r^ + BU,mit)r‘^'‘''^h (shell) (6) 

i ,m 


tj^(r,t) = Aj^(i,m;t)r^Yj^jij(e,<)') , 


*’,m (core) 
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Figure 1. A concentric three fluid 
system used in the theory. 


We now express the boundary functions In the following way: 


(Q*^(t) “ ^ ^ 1 <1 ^ ^ ^ (4(0 

Is} Is} U) ^ ^ ^ 


ii\ u “"ill 

\o) t.m \o| 


.C)Yj^(6,() . 


(11) 


If we substitute these expressions into eqs. (2A) and (3A) , we show that the equations can 
be reduced to an eigenvalue equation which is essentially same as a coupled harmonic 
oscillator; 
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where 
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(15) 

(16) 

(17) 


where x • ^Rq/R^ . Note that the only dimensional quantity is VJ which has the dimension of 
frequency squared. 

From eq. (12), the eigenvalues are given by 




+ 1 , 


so that the normal mode frequencies are given by 


2 K u 
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(18) 


(19) 


From (11) and (17), we obtain the corresponding eigenvectors as 



( 20 ) 


( 21 ) 
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It follows from the orthogonality of eigenvectors 
that 



It is important to note Chat the positive square 
root in ( 20 ) corresponds to the positive square 
root in (18). Consequently for the normal mode 
with the higher frequency, the boundary oscilla- 
tions are in-phase. We call this high frequency, 
mode the "bubble" mode and the lower fre- 
quency, mode the "sloshing" mode. From 

relation ( 22 ), the relative boundary displace- 
ment of the sloshing mode is out of phase. The 
relative boundary displacement of these two modes 
as they were observed in the neutral buoyancy 
tank is shown in Figure 2. The top figure shows 
"bubble" mode in which the two boundaries oscil- 
late in-phase, and the bottom figure shows 
"sloshing" mode in which the two boundaries move 
out of phase. In this compound drop the core and 
the host were silicon oil, and the shell was 
water. 


Velocity potential 


We now obtain expressions of velocity poten- 
tials 1 , 1 , ,s and .0 for the flow f-^Ms in the 
core, shell and the host respectively in terms 
of the interfacial displacements. From eqs. ( 6 ) 
through ( 11 ), the expressions for ,i and 
are obtained as 


( ^ > n, 

1 >m 
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Figure 2. "Bubble" mode (top) and 
"sloshing" mode (bottom) oscilla- 
tions of an oil-water-oil compound 
drop. 


(23) 


where 

s ' and 
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The result for .g, interesting enough, can be expressed as 


(■ s) ' m 




(24) 


(25) 
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where .3 is the shell velocity potential when the core is rigid and vs fs for the case when 

the host is rigid. When exolicitly expressed, they are 


(■ ■ (T-Fttd j<‘ * >>'■ o' * (f) 

(■s)iiB ’ ■ t 1) 0 I*' ^ O I ' (I) + (I) I 


(26) 

( 27 ) 


whore D 


2.+1 .-( 2 :+l) 


(4. 


0 and 


. Of course, conditions 

- 0 were imposed to obtain (26) and (27). 
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Using eqs. (13) Co (19), numerical values of 
the elgenfrequercles for various physical param 
etcrs can be readily obtained. Furthermore, 
simple expressions of several limiting cases as 
shown below can be deduced In a straightforward 
way, 

(I) Simple Drop: We tested the present theory 

by deriving the well-known simple drop case. 

With Ri"0 and following the procedure derived 
in the previous section, we get 

2 . id + !)(■• - 1 ) (> * 2) 

+ ({ + 1 ) . 1 
o 

•> 

which Is Identical to Lamb's result.' 
shows various modes of a simple drop correspond- 
ing to the different values of > . A water drop 
immersed In a neutrally buoyant oil bath was 
axlsymmetrical ly excited by a plunger which was 
In turn connected to a sinusoidally excited loud 
speaker. With l.Scc of drop volume, the inter- 
facial tension, 11.2 dynes/cm, w.is obtained 
using eq. (28) for •' - 2. The higher mode fre- 
quencies obtained using eq . (28) is shown in 
Figure 4 (solid line) and it was compared with 
the actual measurements (circles). The agree- 
ment is excellent. 

(II) Rigid Host ; When the compound drop was 

formed in the rigid host (i.e., , ■ •») , then 

- _ “ 0 and ‘ 


(28) 

Figure 3 





Figure 4. Frequency vs t of a water 
drop in • llicon oil. Drop volume was 
1 . bcc , and » 11.2 dynes /cm as mea- 
sured for i • 2 case. Under these 
conditions, the solid line was 
obtained using eq. (28), and the 
circles were as measured by the 
mechanically excited oscillations. 

-(2t+l), 


(29) 


.-.(c-n .• (( : D(if2)| - , 

Rj |l <c + (<+!) i’j^lt‘’**^+l> + l)(,'-.'j 



If we normalize this with respect to of a 
simple drop which is made out of shell fluid 
immersed in the core fluid, then for 

t-2 and is as shown in Figure 5. Here V'j 

Is the total volume (i.e., Vy-Vshell+Vcore) • 

Nj, rises from zero at Vj/V^.-t, then it 
approaches to Infinity as the core shrinks. As 
the core shrinks with respect to V-j, the core 
motion becomes essentially decoupled from the 
host and approaches to the simple drop fre- 
quency. This can be readily observed from 
eq. (29) as 


(III) Rigid Core ; If the coiiipound drop hac 
rigid core (i.e., .•j’"-), then ._‘0, and 


.'^,(<-1) i (t + l)((+2)h'‘*^ 
Rp{lf.'o+(‘*l). It"'** ♦ «(.> 


.,-(2el), 
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Figure 3. Various modes of a water drop 
oscillations In the silicon oil bath. 

The drop was mechanically excited. 



Of course, when R^-^o, this reduces to eq. (28) for the simple drop. If we let ul this 
simple drop frequency, then the normalized frequency, Nc-u+/uL. for i-2 and pq-p. Is as 
shown in Figure 6. Starting from a very thin shell (i.e., Vt/Vc^ 1), rises sharply to its 
maximum, then gradually approaches to the single drop frequency as Vq'/V(^ increases. The 
agreement with the experimental points is fairly good. 


(iv) Thin Shell : V?hen the shell is very thin (i.e., t ■ 1), then - o, and 


2 ^°o ^ +!)(« +2) 

IiPq + (i + 1) P^l 


(31) 


Note that this is the same expression as that of a simple drop if its Interfacial surface 
tension is replaced by <ro+ •’'i- particular result was also confirmed experimentally in 

the neutral buoyancy tank 



Figure 5. Normalized Frequency, Nq vs 
Vt/Vc when the hose is rigid. 



Figure 6. Normalized frequency, N^, vs 
Vx/Vc when the core is rigid. 


(v) Thick Shell ; When the shell is very 
thick (i.e., !-►"), then we get 

o 4 tto.(i - 1) i (l + l)(t + 2) 

4- [.p ' ^ (V. ly ' pj <32) 

and 

, 4ito„ (i - 1) n a + l)(i + 2) 

— SV 3 (tp-. (IVl) Pi <33) 

where V^. and Vg are the core and the shell 
volumes respectively. These expressions tell 
us that in this limit the "bubble" and 
"sloshing" modes are no longer coupled since 
each mode represents simple drop oscillation. 

(vi) For Pi“Po and oi*oo: When the core and 
the host fluids are the same and the two 
interfacial tensions are same the expressions 
for cannot be made much simpler than 

eq. (19). However, the numerical results of 
these frequencies are shown in Figure 7, 
where Vj-Vc+Vg , and Ns“i*)±/ul where ul is the 
frequency of a simple drop which is made out 
of the same compound drop when the core was 
reduced to zero. In this figure the upper 
five curves represent "bubble" modes for the 
specified values of pi/p when pi>po> and the 
lower three curves represent corresponding 
"sloshing" modes. Experimental points taken 
in a neutral buoyancy tank (pi/p“po/p“D 
show good agreement with the theoretical 
curves. The results for the relative bound- 
ary displacements are shown in Figure 8 for 
an air-water-air type and an oil-water-oil 
type of compound drop systems. We see that 
our preliminary experimental results obtained 
in an neutral buoyancy tank also agree well 
with the theory. 

Remarks on core centering 

An interesting core centering phenomenon 
was observed in our neutral buoyancy experi- 
ments. Initially, a static compound drop 
was prepared so that the inner and the outer 
boundary surfaces were nonconcentric due to 
the slight density mismatch. However, as 
the drop began to oscillate in one of its 
normal mode freauencies (for i-2), the two 
boundaries became concentric within the 
accuracy of oUr observation. Though this 
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Figure 7. Normalized frequency, No, vs Vt/Vq 
when f'i“Co oj[“oo. The upper five curves 
are for the "bubble" modes and the lower 
three curves are for the "Sloshing" modes. 



Figure 8. Relative boundary dis- 
placement of "bubble" modes. Th..' 
upper curve is for an oil-water-oil 
type of compound drop and the lower 
curve is for an air-water-air type 
of compound drop. The experimental 
points are from the neutral buoy- 
ancy experiment. 


centering phenomenon takes place within a few cycles of oscillation, the centering force 
seems to depend, among others, on the oscillation amplitude and the shell thickness. Within 
the approximations used in this paper, the positions of the core anywhere within the drop is 
neutrally stable. It is necessary to consider the next order of approximation to find the 
restoring forces responsible for the centering of the core and the shell. Work on the 
centering phenomenon and its more detailed description vill bi published elsewhere. 


Experimental apparatus and procedure 


(A) Neutral Buoyancy Tank: The heart of the experimental apparatus is a neutral buoy- 

ancy tank which is a lucite box filled with silicon oil (Dow Corning 200, Ics). Taking 
advantage of the fact that the den.sity of silicon oil was less than water, we created a 
vertical density gradient by adding a small amount of freon to the silicon oil until the 
water droplet floated in the middle of the tank. Also imbedded within this box were several 
electrodes for the purpose of excitation and detection of the drop systems. 

(B) Excitation and Detection Procedure: We adopted basically two different detection 

procedures: 1) frequency sweep, and 2) the pulsed transient technique. In the continu- 

ously driving method, the frequency of an oscillating electric field was swept through 
resonances while monitoring the oscillation amplitude using an optical recording and video 
systems (such as fast movie camera or TV video system). This method has the advantage of 
quick identification of the characteristic oscillations and provides detailed information 
about boundary motion. However, it was found that unless the driving field was kept at a 
low level, the resonance frequency measured in this method was e’^roneous. This is due to 
the fact that the large oscillating elect^-ic field changes the surface tension dynamically. 
The detection method which was used extensively in this experiment was the capacitance 
bridge method. The block diagram of this apparatus is shotm in Figure 9. This method is 
similar in principle to more w''''l-known transient techniques in nuclear magnetic resonance 
or optical spectroscopy. The i. sic idea of this technique is to record capacitance varia- 
tion as the drop evolves freely toward equilibrium state from an Initial nonequilibrium 
state. Application of an Intense short electric pulse excites the droplet over a wide 
range of frequency simultaneously. Therefore, the Fourier transformation of this transient 
signal reveals characteristic frequencies simultaneously. 
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Figure 9. A block diagram of our apparatus 
which excites the drop electro-statically 
and detects the ensuing capacitance signal 
caused by the drop deformation. 


Figure 10. Signals obtained by 
the capacitance bridge and their 
Fourier transformed spectra. 


This capacitance bridoe system was constructed arovmd a General Radio Capacitance Bridge 
Unit, a home-build relay system, an A/D converter, and a data processing minicomputer. 

Once the capacitance bridge was balanced with the relay connected to the bridge, an elec- 
tric pulse could be fired by simply controlling the relay switch for the preset time dura- 
tion. Upon the termination of the pulse, the electrodes were switched back to the bridge 
and the ensuing signal was detected and processed. In order to protect the amplifier dur- 
ing the pulse and to shorten the recovery time following the pulse, an amplifier blocking 
relay switch was inserted at the input of the amplifier. 


Figure 10 shows a set of typical data obtained in this method. For the simple drop 
case, the signal is a damped monotonic oscillation which corresponds to a single . esk in 
the frequency domain. Of course, the spectral line width is inversely proportional to the 
damping time constant of the time domain signal. The signal from a comp.yund drop looks 
more complicated. However, Fourier transformation shows two well-resolved peaks each of 
which corresponds to the bubble mode and the slushing mode respectively. Though the capac- 
itance bridge method provides a full spectrum without scanning through resonance frequen- 
cies, it does not provide information about boundary motion. For this reason it is more 
convenient to adopt both continuous and transient methods which will compliment each other. 
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